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We discuss the formation of crystalline electron clusters in semiconductor quantum dots and of 
crystalline patterns of neutral bosons in harmonic traps. In a first example, we use calculations 
for two electrons in an elliptic quantum dot to show that the electrons can localize and form a 
molecular dimer. The calculated singlet-triplet splitting (J) as a function of the magnetic field 
(B) agrees with cotunneling measurements with its behavior reflecting the effective dissociation 
of the dimer for large B. Knowledge of the dot shape and of J{B) allows determination of the 
degree of entanglement. In a second example, we study strongly repelling neutral bosons in two- 
dimensional harmonic traps. Going beyond the Gross-Pitaevskii (GP) mean-field approximation, we 
show that bosons can localize and form polygonal-ring-like crystalline patterns. The total energy 
of the crystalline phase saturates in contrast to the GP solution, and its spatial extent becomes 
smaller than that of the GP condensate. 
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Explorations of the size-dependent evolution of the 
properties of materials are at the frontier of modern con- 
densed matter and materials research. Indeed, investiga- 
tions of clusters containing a finite, well-defined, number 
of elementary building units (atoms, molecules, electrons 
or other elementary constituents), allow investigations 
of the transition from the atomic or molecular regime 
to the finite nano-aggregate domain, and ultimately of 
the convergence with increasing size to the condensed 
phase, extended system, category. Moreover, investiga- 
tions of clusters provide opportunities for discovery of 
novel properties and phenomena that are intrinsic prop- 
erties of finite systems, distinguishing them from bulk 
materials Q. 

Commonly, studies of materials clusters involve 
atoms or molecules interacting through electrostatic or 
electromagnetic potentials, with the heavier nuclei being 
the "structural skeleton" and the much lighter electrons 
serving as the "glue" that binds the atoms together. In 
this article we focus on novel, somewhat exotic, types of 
clusters. In particular, we discuss clusters of electrons 
in man-made (artificial) quantum dots (QDs) created 
through lithographic and gate-voltage techniques at 
semiconductor interfaces, and clusters of neutral atoms 
in traps under conditions that may relate to formation 
of Bose-Einstein condensates (BECs). We illustrate 
that these cluster systems reveal interesting emergent 
physical behavior arising from spontaneous breaking of 
spatial symmetries; symmetry breaking (SB) is defined 
as a circumstance where a lower energy solution of the 
Schrodinger equation is found that is characterized by 
a lower symmetry than that of the Hamiltonian of the 
system. Such SB is exhibited through the formation 
of clusters of localized electrons (often called Wigner 
molecules, WMs) in two-dimensional (2D) QDs (see Fig. 
1 and in particular our discussion of two-electron WMs 
below 0). Symmetry breaking is also manifested in 



the transition 0, induced by increasing the interatomic 
repulsive interaction strength, of the BEC state of 
neutral atoms confined by a parabolic 2D trap to a 
crystalline cluster state. 



Two-electron Wigner molecules. Electron localiza- 
tion leading to formation of molecular-like structures 
[the aforementioned Wigner molecules] within a sin- 
gle circular two-dimensional quantum dot at zero mag- 
netic field (B) has been theoretically predicted to occur 
H H i, Q, H H H EH, as the strength of the e-e re- 
pulsive interaction relative to the zero-point energy in- 
creases, as expressed through an inceasing value of the 




FIG. 1: Unrestricted Hartree-Fock electron density in a 
parabolic QD for TV = 19 electrons and S z = 19/2, exhibit- 
ing breaking of the circular symmetry at Rw = 5 and zero 
magnetic field. Remaining parameters are: parabolic confine- 
ment, %u)o = 5 meV; effective mass m* = 0.067m e . Distances 
are in nanometers and the electron density in 10 -4 nm -2 . 
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Wigner parameter Rw, defined as Rw = Z 2 e 2 / '(huiolo) , 
with Iq — \/~Kj (mwo) being the characteristic harmonic- 
oscillator length of the confining potential. (The sub- 
script W in the case of a Coulomb force stands for 
"Wigner", since the confined clusters of localized elec- 
trons may be viewed as finite-size precursors of the bulk 
Wigner crystal H3-) Formation of such "molecular struc- 
tures" is a manifestation of spontaneous symmetry break- 
ing associated with a quantum phase transition, occuring 
at zero magnetic field for circular 2D QDs for Rw > 1 
and involving the crossover from a liquid-like state to a 
crystalline one. This crossover, described using 
try breaking [J] , was confirmed in other studies 
using a variety of methods. The degree of electron lo- 
calization, that underlies the appearance of crystalline 
patterns, has been described |jj as a progression from 
"weak" to "strong" Wigner molecules as a function of in- 
creasing Rw (or equivalently decreasing density). In high 
magnetic fields, the rotat ing electron molecules exhibit 
magic angular momenta [Til IbH fisl flfij ] corresponding 
to fractional quantum Hall effect (FQHE) fillings. This 
has led [lj, Il5l |l6( to the derivation of an analytic trial 
wave function that provides a better description of the 
finite-size analogs of the FQHE in comparison with the 
Laughlin and composite-fermion |l8j | wave functions. 

Because of the finite size, Wigner molecules are ex- 
pected to show new behavior that differs from the clas- 
sical Wigner crystal familiar from Solid State Physics. 
The limit of a classical Wigner crystal is expected to 
be reached for a higher number of electrons N and very 
large Rw 4J. In the following we use the term "Wigner 
molecule" even in the case of only two localized electrons. 
For this case the WM exhibits close analogies to an H2 
natural molecule, as described below. 

Here we focus on a two-electron (2e) WM, in light of 
the current experimental effort [h!.[2o| aiming at imple- 
mentation of a spin-based [ilj solid-state quantum logic 
gate that employs two coupled one-electron QDs (dou- 
ble dot). We present an exact diagonalization (EXD) 
and an approximate (generalized Heitler-London, GHL) 
microscopic treatment for two electrons in a single el- 
liptic QD specified by the parameters of a recently in- 
vestigated experimental device jl^. While formation of 
Wigner molecules in circular QDs requires weak confine- 
ment (that is small loq in the expression for Rw given 
above), and thus large dots of lower densities (so that 
the interelectron repulsion dominates) , we show that for- 
mation of such WMs is markedly enhanced in highly 
deformed (e.g., elliptic) dots due to their lower symme- 
try. The calculations provide a good description of the 
measured J(B) curve (the singlet-triplet splitting) when 
screening [23|, |2J] due to the metal gates and leads is in- 
cluded (in addition to the dielectric constant of the semi- 
conductor, GaAs). In particular, our results reproduce 
the salient experimental findings pertaining to the van- 
ishing of J(B) for a finite value of B ~ 1.3 T [associated 



with a change in sign of J(B) indicating a singlet-triplet 
(ST) transition], as well as the flattening of the J(B) 
curve after the ST crossing. These properties, and in par- 
ticular the latter one, are related directly to the formation 
of an electron molecular dimer and its effective dissocia- 
tion for large magnetic fields. The effective dissociation 
of the electron dimer is most naturally described through 
the GHL approximation, and it is fully supported by the 
more accurate, but physically less transparent, EXD. 

Of special interest for quantum computing [2l| is the 
degree of entanglement exhibited by the two-electron 
molecule in its singlet state. Entanglement is a purely 
quantum mechanical phenomenon in which the quantum 
state of two or more objects cannot be described indepen- 
dently of each other, even when the individual objects are 
spatially separated. The highest degree of entanglement 
occurs at full separation, as discussed in the celebrated 
EPR (Einstein, Podolsky, and Rosen) paper [ilj . Elec- 
trons confined in a quantum dot are not necessarily spa- 
tially separated from each other, and consequently their 
degree of entanglement may be lower than the maximal 
one, as shown by us below. 

Here, in relation to the microscopic calculations, we 
investigate two different measures of entanglement. The 
first, known as the concurrence (C) for two indistinguish- 
able fermions [2H. , has been used in the analysis of 
the experiment in Ref. (this measure is related to 
the operational cycle of a two-spin-qubit quantum logic 
gate H^, The second measure, referred to as the 

von Neumann entropy (S) for indistinguishable particles, 
has been developed in Ref. |2g and used in Ref. [29| . 
We show that the present wave-function-based methods, 
in conjunction with the knowledge of the dot shape and 
the J(B) curve, enable theoretical determination of the 
degree of entanglement, in particular for the elliptic QD 
of Ref. [2^ ■ The increase in the degree of entanglement 
(for both measures) with stronger magnetic fields 
correlates with the dissociation of the 2e molecule. This 
supports the experimental assertion 22] that cotunncl- 
ing spectroscopy can probe properties of the electronic 
wave function of the QD, and not merely its low-energy 
spectrum. Our methodology can be straightforwardly 
applied to other cases of strongly-interacting devices, 
e.g., double dots with strong interdot-tunneling. 



Clusters of neutral bosons in harmonic traps. 

Bose-Einstein condensates (BECs) in harmonic traps 
[iffll IsH ] are normally associated with weakly interact- 
ing neutral atoms, and their physics is described ade- 
quately by the Gross-Pitaevskii (GP) mean- field theory 
|32j. Lately, however, experimental advances in control- 
ling the interaction strength [sjj 0, I^H |^] permit the 
production of novel bosonic states in the regime of strong 
interparticle repulsions. Theoretical efforts motivated 
by this capability include studies of the Bose-Hubbard 
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model |3jj, |38|], and investigations about the "fermion- 
ization" limit of an one-dimensional (ID) gas of trapped 
impenetrable bosons [H 0, lili. o ften referred to as the 
Tonks-Girardeau (TG) regime [3^. E3 ] . Here we address 
the problem of strongly repelling (impenetrable) bosons 
in higher dimensions. In particular, we discuss 2D in- 
teracting bosons in a circular harmonic trap, with the 
extension to 3D systems being straightforward. To this 
end, we use computational methods that go beyond the 
GP method. 

We explore the transition from a BEC (diffuse cloud) 
state to a crystalline phase, in which the trapped local- 
ized bosons form crystalline patterns. At the mean-field 
level, these crystallites are static and are portrayed di- 
rectly in the single-particle densities. After restoration 
of rotational symmetry, the single-particle densities are 
circularly symmetric, and thus the crystalline symmetry 
becomes "hidden" ; however, it can be revealed in the con- 
ditional probability distribution (CPD, anisotropic pair 
correlation), P(r, ro), which expresses the probability of 
finding a particle at r given that the "observer" (i.e. , ref- 
erence point) is riding on another particle at r [lCillfij. 

METHODS 

Two-electron quantum dot: Microscopic treat- 
ment. The Hamiltonian for two 2D interacting electrons 
is 

H = H(r 1 )+H(r 2 )+e 2 /( K r 12 ), (1) 

where the last term is the Coulomb repulsion, k is the 
dielectric constant, and 7-12 = ri — r2 1 . H(r) is the 
single-particle Hamiltonian for an electron in an exter- 
nal perpendicular magnetic field B and an appropriate 
confinement potential. When position-dependent screen- 
ing is included, the last term in Eq. is modified by 
a function of r\ 2 (see below). For an elliptic QD, the 
single-particle Hamiltonian is written as 

H(v) =T+ 1 -m* + coy) + ^B • s, (2) 

where T = (p - eA/c) 2 /2m* , with A = 0.5(-By, Bx, 0) 
being the vector potential in the symmetric gauge, m* 
is the effective mass and p is the linear momentum of 
the electron. The second term is the external confining 
potential; the last term is the Zeeman interaction with 
g* being the effective g factor, [Ib the Bohr magneton, 
and s the spin of an individual electron. 

The GHL method for solving the Hamiltoninian (JTJ 
consists of two steps. In the first step, we solve selfcon- 
sistcntly the ensuing unrestricted Hartree-Fock (UHF) 
equations allowing for lifting of the double-occupancy re- 
quirement (imposing this requirement gives the restricted 
HF method, RHF). For the S z = solution, this step 



produces two single-electron orbitals ul,r{y) that are lo- 
calized left (L) and right (R) of the center of the QD [un- 
like the RHF method that gives a single doubly-occupied 
elliptic (and symmetric about the origin) orbital]. At 
this step, the many-body wave function is a single Slater 
determinant \E'uhf(1 t,2 J.) = \ul(1 *])u r (2 I)) made 
out of the two occupied UHF spin-orbitals f) = 

Ui(ri)a(l) and u R {2 [) = u R (r 2 )/3(2), where a(0) de- 
notes the up (down) [\ (J,)] spin. This UHF determinant 
is an eigenfunction of the projection S z of the total spin 
S = si + S2, but not of S 2 (or the parity space-reflection 
operator) . 

In the second step, we restore the broken parity and 
total-spin symmetries by appl ying to the UHF determi- 
nant the projection operator |43l E3] P s,t = 1 =p tui2, 
where the operator zvi 2 interchanges the spins of the 
two electrons; the upper (minus) sign corresponds to the 
singlet. The final result is a generalized Heitler-London 
(GHL) two-electron wave function ^qhlC 1 "!' r 2) f° r tne 
ground-state singlet (index s) and first-excited triplet (in- 
dex t), which uses the UHF localized orbitals, 

*^ L (r!,r 2 ) oc (u L (r 1 )u R {r 2 )±u L (r 2 )u R (r 1 ))x s '\ (3) 

where x s,t — ( a (l)/?(2) T ck(2)/?(1) ) is the spin function 
for the 2e singlet and triplet states. The general for- 
malism of the 2D UHF equations and of the subsequent 
restoration of broken spin symmetries can be found in 

Refs. Ill El BE!. 

The use of optimized UHF orbitals in the GHL is suit- 
able for treating single elongated QDs. The GHL is 
equally applicable to double QDs with arbitrary interdot- 
tunneling coupling 0, In contrast, the Heitler- 

London (HL) treatment [46j (known also as Valence 
bond) , where non-optimized "atomic" orbitals of two iso- 
lated QDs are used, is appropriate only for the case of a 
double dot with small interdot-tunneling coupling 21]. 

The orbitals ul, R (y) are expanded in a real Cartesian 
harmonic-oscillator basis, i.e., 

K 

u L , R (r) = J2 C f' R Vj(r), ( 4 ) 

where the index j = (m,n) and <pj(r) = X m [x)Y n {y), 
with X m (Y n ) being the eigenfunctions of the one- 
dimensional oscillator in the x(y) direction with fre- 
quency cu x (ujy). The parity operator V yields VX m (x) — 
(—l) m X m (x), and similarly for Y n (y). The expansion co- 
efficients Cj ' are real for B = and complex for finite 
B. In the calculations we use K = 79, yielding conver- 
gent results. 

In the exact-diagonalization method, the many-body 
wave function is written as a linear superposition over 
the basis of non-interacting two-electron determinants, 
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i.e., 

2K 

*&D(ri,rs) =X;^/lV-(l;^(2;i)), (5) 

where ip{l;i) = <fi(l f) if 1 < i < K and 
■0(1; i) = <Pi-KO- I) if K + 1 < % < 2K [and simi- 
larly for 0(2, j)]. The total energies -Eg XD an< ^ * ne 
coefficients fiy are obtained through a "brute force" 
diagonalization of the matrix eigenvalue equation corre- 
sponding to the Hamiltonian in Eq. QJ. The EXD wave 
function does not immediately reveal any particular 
form, although, our calculations below show that it can 
be approximated by a GHL wave function in the case of 
the elliptic dot under consideration. 

Two-electron quantum dot: Measures of entan- 
glement. To calculate the concurrence C [2H. |27|. one 
needs a decomposition of the GHL wave function into 
a linear superposition of orthogonal Slater determinants. 
Thus one needs to expand the nonorthogonal u L,R (r) or- 
bitals as a superposition of two other orthogonal ones. 
To this effect, we write u LR (r) ex $ + (r) ±£$~(r), where 
$ + (r) and 3> - (r) are the parity symmetric and antisym- 
metric (along the ir-axis) components in their expansion 
given by Eq. Q. Subsequently, with the use of Eq. J3J), 
the GHL singlet can be rearranged as follows: 

*&hl « |$+(1 T)$+(2 I)) - t)*-(2 I)), (6) 

where the so-called interaction parameter [23 , rj — £ 2 , is 
the coefficient in front of the second determinant. Know- 
ing X] allows a direct evaluation of the concurrence of the 
singlet state, since C s = 2-q/(l+jf) [13 . Note that $+(r) 
and $ _ (r) are properly normalized. It is straightforward 
to show that 77 = (1 - \S LR \)/(1 + \S L r\), where S LR 
(with \Slr\ < 1) is the overlap of the original u L,R (r) 
orbitals. 

For the GHL triplet, one obtains an expression inde- 
pendent of the interaction parameter 77, i.e., 

*ghl « |*+(1 T)*~(2 I)) + |$ + (1 I)$-(2 |)>, (7) 

which is a maximally (C* = 1) entangled state. Note that 
underlying the analysis of the experiments in Ref. 22] is 
a conjecture that wave functions of the form given in Eqs. 
© and @ describe the two electrons in the elliptic QD. 

To compute the von Neumann entropy, one needs to 
bring both the EXD and the GHL wave functions into 
a diagonal form (the socalled "canonical form" [2^, ) > 
i.e., 

M 

*^ D (n,r 2 ) = 2^ - 1)*(2; 2k)), (8) 

k=l 

with the $(i)'s being appropriate spin orbitals resulting 
from a unitary transformation of the basis spin orbitals 



0(j)'s [see Eq. J5J]; only terms with Zfc / contribute. 
The upper bound M can be smaller (but not larger) than 
K (the dimension of the single-particle basis); M is re- 
ferred to as the Slater rank. One obtains the coefficients 
of the canonical expansion from the fact that the \z^\ 2 are 
eigenvalues of the hermitian matrix 17+0 [fi, see Eq. |JS}, 
is antisymmetric] . The von Neumann entropy is given by 

M 

S = -5> fc | 2 log 2 (|z fe | 2 ) (9) 

k=l 

with the normalization J\=i \ z k\ 2 = 1< Note that 
the GHL wave functions in Eqs. © and J7J are 
already in canonical form, which shows that they 
always have a Slater rank of M = 2. One finds 
5 GHL = l°g 2 (l + V 2 ) - T7 2 log 2 (77 2 )/(l + V 2 ), and 
<Sq HL = 1 for all B. For large B, the overlap between 
the two electrons of the dissociated dimer vanishes, and 
thus 77 — ► 1 and 5q HL — > 1. 

Neutral repelling bosons in harmonic traps: 
Gross-Pitaevskii equation. Mean-field symmetry 
breaking for bosonic systems has been discussed earlier 
in the context of two-component condensates, where each 
species is associated with a different space orbital |48u49| . 
We consider here one species of bosons, but allow each 
particle to occupy a different space orbital 4>i(vi). The 
permanent |$jv) = Perm [0i(ri ),..., 4>n(^n)\ serves as 
the many-body wave function of the unrestricted Bosc- 
Hartree-Fock (UBHF) approximation. This wave func- 
tion reduces to the Gross-Pitaevskii form with the re- 
striction that all bosons occupy the same orbital 4>q(t), 

i.e., |$^ p ) = rifc=i0o( r i)j an d 0o( r ) is determined self- 
consistently at the restricted Bose-Hartree-Fock (RBHF) 
level via the equation |5(| 

[flo(ri) + {N-l) J dr 2 05(r 2 )V(ri, r 2 )0 o (r 2 )]0o(ri) 
= eo0o(ri). (10) 

Here U(ri,r 2 ) is the two-body repulsive interaction, 
which is taken to be a contact potential, Vs — gS(ri — r 2 ), 
for neutral bosons. The single-particle hamiltonian is 
given by Hq(t) = — h 2 V 2 / (2m) + muj 2 r 2 /2, where ojq 
characterizes the harmonic confinement. 



Neutral repelling bosons in harmonic traps: Sym- 
metry breaking. We simplify the solution of the UBHF 
problem by considering explicit analytic expressions for 
the space orbitals 4>i(vi). In particular, since the bosons 
must avoid occupying the same position in space in order 
to minimize their mutual repulsion, we take all the or- 
bitals to be of the form of displaced Gaussians, namely, 
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(f>i{vi) = n~ 1 / 2 <7~ 1 exp[— (rj - a 4 ) 2 /(2(7 2 )]. The positions 
aj describe the vertices of concentric regular polygons, 
with both the width a and the radius a = |aj | of the regu- 
lar polygons determined variationally through minimiza- 
tion of the total energy -Eubhf = (^n\H \$n) /(&n\$n), 

where H = Yn=i H o( r i) + T,f<j v ( T i> T j) is the many- 
body hamiltonian. 

With the above choice of localized orbitals, the 
unrestricted permanent |$jv) breaks the continuous 
rotational symmetry. However, the resulting energy gain 
becomes substantial for stronger repulsion. Controlling 
this energy gain (the strength of correlations) is the 
ratio Rs between the strength of the repulsive potential 
and the zero-point kinetic energy. Specifically, for a 2D 
trap, one has Rs = gm/ (2irh 2 ) for a contact potential. 

Neutral repelling bosons in harmonic traps: 
Restoration of broken symmetry. Although the op- 
timized UBHF permanent \$n) performs exceptionally 
well regarding the total energies of the trapped bosons, 
in particular in comparison to the resctricted wave func- 
tions (e.g., the GP anzatz), it is still incomplete. In- 
deed, due to its localized orbitals, |$jv) does not preserve 
the circular (rotational) symmetry of the 2D many-body 
hamiltonian H . Instead, it exhibits a lower point-group 
symmetry, i.e., a C2 symmetry for N = 2 and a C5 one 
for the (1,5) structure of N = 6 (see below). As a re- 
sult, |$jv) does not have a good total angular momen- 
tum. This paradox is resolved through a post-Hartree- 
Fock step of restoration of broken symmetries via pro- 
jection techniques [isl l5ll|. yielding a new wave function 
l^jvx) w ith a definite angular momentum L, that is 

2n\^l)= T d^i^e^, (11) 
Jo 

where \&N("f)} is the original UBHF permanent having 
each localized orbital rotated by an azimuthal angle 7, 
with L being the total angular momentum. The projec- 
tion yields wave functions for a whole rotational band. 
Note that the projected wave function in Eq. l(TT|) 

may be regarded as a superposition of the rotated per- 
manents \(&n ("/)}, thus corresponding to a "continuous- 
configuration- interaction" solution. 

Here, we are interested in the projected ground-state 
(L = 0) energy, which is given by 

< RJ = (<y\m^)/(*Z J \^)- (12) 

RESULTS AND DISCUSSION 

Two-electron quantum dot. To model the experi- 
mental elliptic QD device, we take, following Ref. [22| . 
Tilj x = 1.2 mcV and fnj v — 3.3 meV. The effective mass 
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FIG. 2: The singlet-triplet splitting J = E s - E l as a func- 
tion of the magnetic field B for an elliptic QD with hui x = 1.2 
meV and hui y = 3.3 meV (these values correspond to the 
device of Ref. H2)- Solid line: GHL (broken-symmetry 
UHF + restoration of symmetries) results with a coordinate- 
independent screening (re = 22). Dashed line: EXD results 
with re = 12.9 (GaAs), but including screening with a coor- 
dinate dependence according to Ref. |24j| and d = 18.0 nm 
(see text). The rest of the material parameters used are: 
m*(GaAs)= 0.067m e , and g* = (see text). The experi- 
mental measurements |22| are denoted by open squares. Our 
sign convention for J is opposite to that in Ref. |22|. 



of the electron is taken as m* — 0.067m e (GaAs). Since 
the experiment did not resolve the lifting of the triplet 
degeneracy caused by the Zeeman term, we take g* = 0. 
Using the two step method, we calculate the GHL singlet- 
triplet splitting Jqhl(-B) = E s wh (B) - E^^B) as a 
function of the magnetic field in the range < B < 2.5 
T. Screening of the e-e interaction due to the metal gates 
and leads must be considered in order to reproduce the 
experimental J(B) curve. This screening can be mod- 
eled, to first approximation, by a position-independent 
adjustment of the dielectric constant k j5^. Indeed, with 
k = 22.0 (instead of the GaAs dielectric constant, i.e., 
k = 12.9), good agreement with the experimental data is 
obtained [see Fig. 2]. In particular, we note the singlet- 
triplet crossing for B » 1.3 T, and the flattening of the 
J(B) curve beyond this crossing. 

We have also explored, particularly in the context of 
the EXD treatment, a position-dependent screening us- 
ing the functional form, (e 2 /Kri 2 )[l — (1 + 4d 2 /r 2 2 ) _1//2 ], 
proposed in Ref. |24|. with d as a fitting parameter. The 
Jexd(B) result for d = 18.0 nm is depicted in Fig. 2 
(dotted line), and it is in very good agreement with the 
experimental measurement. 

The singlet state electron densities from the GHL and 
the EXD treatments at B = and B = 2.5 T are dis- 
played in Fig. 3. These densities illustrate the dissocia- 
tion of the electron dimer with increasing magnetic field. 
The asymptotic convergence (beyond the ST point) of the 
energies of the singlet and triplet states, i.e., [J(B) — > 
as B — > 00] is a reflection of the dissociation of the 2e 
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(a) 



GHL B=0 




(b) 



EXD B=0 




FIG. 3: Total electron densities (EDs) associated with the 
singlet state of the elliptic dot at B — and B = 2.5 T. (a) 
The GHL densities, (b) The EXD densities. The rest of the 
parameters and the screening of the Coulomb interaction are 
as in Fig. 2. Lengths in nm and densities in 10 -4 nm~ 2 . 



molecule, since the ground-state energy of two fully spa- 
tially separated electrons (zero overlap) does not depend 
on the total spin. 

In contrast, the singlet-state RHF electron densities 
fail to exhibit formation of an electron dimer for all values 
of B. This underlies the failure of the RHF method to 
describe the behavior of the experimental J(B) curve. 
In particular, Jr,hf(-E> = 0) has the wrong sign, while 
Jrhf(B) diverges for high B as is the case for the RHF 
treatment of double dots (see Ref. |43| ). 

For the GHL singlet, using the overlaps of the left 
and right orbitals, we find that starting with rj = 0.46 
(C s = 0.76) at B = 0, the interaction parameter (singlet- 
state concurrence) increases monotonically to ry = 0.65 
(C s = 0.92) at B = 2.5 T. At the intermediate value 
corresponding to the ST transition (B = 1.3 T), we find 
r\ = 0.54 (C s = 0.83). 

Our B = theoretical results for rj and C s are in re- 
markable agreement with the experimental estimates [22j | 
of i] = 0.5 ± 0.1 and C s = 0.8, which were based solely 
on conductance measurements below the ST transition 
(i.e., near B — 0). We note that, for the RHF, Cf> HF = 0, 
since a single determinant is unentangled for both the 
two measures considered here. 

Since the EXD singlet has obviously a Slater rank 
M > 2, the definition of concurrence is not applicable 
to it. The von Neumann entropy for the EXD singlet 
(<5g XD ) is displayed in Fig. 4, along with that (<Sq HL ) 
of the GHL singlet. S^ XD and Sq HL are rather close to 
each other for the entire B range, and it is remarkable 
that both remain close to unity for large B, although 
the maximum allowed mathematical value is \og 2 (K) 
[as aforementioned we use K = 79, i.e., log 2 (79) = 6.3]; 
this maximal value applies for both the EXD and GHL 
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FIG. 4: Von Neumann entropy for the singlet state of the 
elliptic dot as a function of the magnetic field B. Solid line: 
GHL. Dashed line: EXD. The rest of the parameters and the 
screening of the Coulomb interaction are as in Fig. 2. At the 
top, we show histograms for the |zfc| 2 coefficients [sec Eq. JjJ] 
of the singlet state at B — 1.3 T, illustrating the dominance 
of two configurations. Note the small third coefficient |z3| 2 = 
0.023 in the EXD case. 



approaches. The saturation of the entropy for large 
B to a value close to unity reflects the dominant (and 
roughly equal at large B) weight of two configurations 
in the canonical expansion [see Eq. ©] of the EXD 
wave function, which are related to the two terms 
(M = 2) in the canonical expansion of the GHL singlet 
[Eq. This is illustrated by the histograms of the 

|z|| 2 coefficients for B = 1.3 T at the top of Fig. 4. 
These observations support the GHL approximation, 
which is computationally less demanding than the exact 
diagonalization, and can be used easily for larger TV. 



Neutral repelling bosons in harmonic traps. In 

Fig. 5, we display as a function of the parameters Rg 
the total energies for TV = 6 bosons calculated at several 
levels of approximation. In both cases the lowest UBHF 
energies correspond to a (1,5) crystalline configuration, 
namely one boson is at the center and the rest form a 
regular pentagon of radius a. Observe that the GP total 
energies are slightly lower than the ^rbhf ones ; however, 
both exhibit an unphysical behavior since they diverge as 
Rs — ► oo. This behavior contrasts sharply with that of 
the unrestricted Hartree-Fock energies, -Etjbhf and PRJ 
(see below), which saturate as Rs — ► oo; in fact, a value 
close to saturation is achieved already for Rs ~ 10. We 
have checked that for all cases with N — 2 — 7, the to- 
tal energies exhibit a similar behavior. For a repulsive 




contact 
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UBHF 

" PRJ 



100 _, 200 



FIG. 5: Total energies as a function of Rg for various approx- 
imation levels, calculated for N = 6 harmonically confined 
2D bosons in the (1,5) lowest-energy configuration. Nota- 
tion: RBHF/G - Restricted Bose-Hartree-Fock (RBHF) en- 
ergy, Srbhf , with the common orbital 0o(r) approximated 
by a Gaussian centered at the trap origin; GP - the Gross- 
Pitaevskii energy; PRJ - the energy of the symmetry-restored 
state obtained via projection of the (unrestricted) UBHF 
state. Energies in units of hu>o- 



contact potential, the saturation of the UBHF energies 
is associated with the ability of the trapped bosons (in- 
dependent of N) to minimize their mutual repulsion by 
occupying different positions in space, and this is one of 
our central results. For N — 2, the two bosons localize 
at a distance 2a apart to form an antipodal dimer. For 
N < 5 the preferred UBHF crystalline arrangement is a 
single ring with no boson at the center [usually denoted 
as (0, N)]. N = 6 is the first case having one boson 
at the center [designated as (l,N — 1)], and the (0,6) 
arrangement is a higher energy isomer. The structural 
parameters (e.g., the width of the Gaussian orbitals and 
the radii of the polygonal ring, calculated via the UBHF 
method, show a saturation behavior similar to that illus- 
trated above for the energy of the system |3(. In contrast, 
the width of the condensate cloud (i.e., the GP solution) 
diverges with increasing repulsion strength (R$). 

The saturation found here for 2D trapped bosons inter- 
acting through strong repelling contact potentials is an 
illustration of the "fermionization" analogies that appear 
in strongly correlated systems in all three dimensionali- 
ties. Indeed such e nerg y saturation has been shown for 
the TG ID gas [3^, l42| . and has also been discussed for 
certain 3D systems ( i.e., three trapped bosons and 
an infinite boson gas |54|). Saturation of the energy and 
the length of the trapped atom cloud (and thus of the in- 
terparticle distance) has been measured recently for the 
ID TG gas (see in particular Fig. 3 and Fig. 4 in Ref. 
|36| and compare to the similar trends predicted here for 
the 2D case in Fig. 5). 

For N = 6 2D bosons, Fig. 5 shows that the E% RJ 
energies share with the UBHF ones the saturation prop- 
erty for the case of a contact-potential repulsion. How- 
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FIG. 6: (a-c): Single-particle densities for N = 6 2D har- 
monically trapped neutral bosons with a contact interaction 
and Rg — 25. (a) The single-orbital self-consistent GP case, 
(b) The symmetry broken UBHF case (static crystallite), (c) 
The projected case [symmetry-restored wave function, see Eq. 
Illl l. The crystalline structure of the outer ring in this last 
case is "hidden", but it is revealed in the conditional proba- 
bility distribution (ToL ITr| displayed in (d) , where the obser- 
vation point is denoted by a black dot (on the right). Lengths 
in units of Iq. 

ever, the projection brings further lowering of the total 
energies compared to the UBHF ones. (The projected 
ground state is always lower in energy than the original 
broken-symmetry one |5o|.) Thus, for strong interactions 
(large values of R$) the restoration-of- broken-symmetry 
step yields an excellent approximation of both the exact 
many-body wave function and the exact total energy. 

The transformations of the single-particle densities 
(displayed in Fig. 6 for N = 6 neutral bosons inter- 
acting via a contact potential and Rs — 25) obtained 
from application of the successive approximations pro- 
vide an illustration of the two-step method of symme- 
try breaking with subsequent symmetry restoration. In- 
deed, the GP single-particle density [Fig. 6(a)] is circu- 
larly symmetric, but the UBHF one [Fig. 6(b)] explic- 
itly exhibits a (1,5) crystalline configuration. After sym- 
metry restoration [Fig. 6(c)], the circular symmetry is 
re-established, but the single-particle density is radially 
modulated unlike the GP density. In addition, the crys- 
talline structure in the projected wave function is now 
hidden; however, it can be revealed through the use of the 
CPD 10, 16| [see Fig. 6(d)], which resembles the (crys- 
talline) UBHF single-particle density, but with one of 
the humps on the outer ring missing (where the observer 
is located). In particular, P(ro,ro) « and the boson 
associated with the observer is surrounded by a "hole" 
similar to the exchange-correlation hole in electronic sys- 
tems. This is another manifestation of the "fermion- 
ization" of the strongly repelling 2D bosons. However, 
here as in the ID TG case |39|, |42J, the vanishing of 
P(ro,ro) results from the impenetrability of the bosons. 
For the GP condensate, the CPD is independent of ro, 
i.e., PGp( r i r o) l^oMI 2 ; reflecting the absence of any 
space correlations. 

It is of importance to observe that the radius of the 
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BEC [GP case, Fig. 6(a)] is significantly larger than the 
actual radius of the strongly-interacting crystalline phase 
[projected wave function, Fig. 6(c)]. This is because the 
extent of the crystalline phase saturates, while that of 
the GP condensate grows with no bounds as Rg — » oo. 
Such dissimilarity in size (between the condensate and 
the strongly-interacting phase) has been also predicted 
[ioj l for the trapped ID Tonks- Girardeau gas and indeed 
observed experimentally j^] • In addition, the 2D single- 
particle momentum distributions for neutral bosons have 
a one-hump shape with a maximum at the origin (a be- 
havior exhibited also by the trapped ID TG gas). The 
width of these momentum distributions versus R$ in- 
creases and saturates to a finite value, while that of the 
GP solution vanishes as R$ — > oo. 



SUMMARY 

In this paper, we explored symmetry-breaking transi- 
tions predicted to occur in confined fcrmionic and bosonic 
systems when the strength of the interparticle repulsive 
interactions exceeds an energy scale that characterizes 
the degree of confinement. For two electrons in an el- 
liptic QD, we predicted formation and effective dissocia- 
tion (with increasing magnetic field) of an electron dimer, 
which is reflected in the behavior of the computed singlet- 
triplet splitting, J(B), that agrees well (Fig. 2) with mea- 
surements |22j. 

Furthermore, we showed that, from a knowledge of the 
dot shape and of J(B), theoretical analysis along the lines 
introduced here allows probing of the correlated ground- 
state wave function and determination of its degree of 
entanglement. This presents an alternative to the exper- 
imental study where determination of the concurrence 
utilized conductance data |22j . Such information is of 
interest to the implementation of spin-based solid-state 
quantum logic gates. 

For the case of 2D trapped bosonic clusters, we found 
with increasing repulsive two-body interaction localiza- 
tion of the bosons in the trap, resulting in formation of 
crystalline patterns made of polygonal rings; while we 
have focused here on repulsive contact interactions, sim- 
ilar results were obtained also for a Coulomb repulsion 
0. 

These results provide the impetus for experimental ef- 
forts to access the regime of strongly repelling neutral 
bosons in two dimensions. To this end we anticipate that 
extensions of methodologies developed for the recent re- 
alization of the Tonks-Girardeau regime in ID (using a 
finite small number of trapped 87 Rb and optical lattices, 
with a demonstrated wide variation of Rs from 5 to 200 
and from 1 to 5 :36]) will prove most promising. Con- 
trol of the interaction strength via the use of the Feshbach 
resonance may also be considered |33j |. 

This research is supported by the US D.O.E. (Grant 
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